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Abstract. We present 325 MHz (90 cm wavelength) radio observations of 
ultracool dwarfs TVLM 513-46546 and 2MASS J0036+1821104 using the Very 
Large Array (VLA) in June 2007. Ultracool dwarfs are expected to be 
undetectable at radio frequencies, yet observations at 8.5 GHz (3.5 cm) and 
4.9 GHz (6 cm) of have revealed sources with > 100 /ijy quiescent radio flux 
and > 1 mjy pulses coincident with stellar rotation. The anomalous emission 
is likely a combination of gyrosynchrotron and cyclotron maser processes in 
a long-duration, large-scale magnetic field. Since the characteristic frequency 
for each process scales directly with the magnetic field magnitude, emission at 
lower frequencies may be detectable from regions with weaker field strength. 
We detect no significant radio emission at 325 MHz from TVLM 513-46546 or 
2MASS J0036+1821104 over multiple stellar rotations, establishing 2.5ct total 
flux limits of 795 fi3y and 942 fj,Jj respectively. Analysis of an archival VLA 
1.4 GHz observation of 2MASS J0036+1821104 from January 2005 also yields a 
non-detection at the level of < 130 /^Jy . The combined radio observation history 
(0.3 GHz to 8.5 GHz) for these sources suggests a continuum emission spectrum 
for ultracool dwarfs which is either flat or inverted below 2-3 GHz. Further, if the 
cyclotron maser instability is responsible for the pulsed radio emission observed on 
some ultracool dwarfs, our low-frequency non-detections suggest that the active 
region responsible for the high-frequency bursts is confined within 2 stellar radii 
and driven by electron beams with energies less than 5 keV. 
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1. Introduction 



Ultracool dwarfs (UCDs) describe a subsection of stellar objects located on the 
boundary between more massive stars and sub-stellar bodies such as gas giant planets. 
It includes fully convective, very low mass stars (M7 and lower) and all brown dwarfs. 
X-ray and Ha intensity (typical proxies for ma gnetic activity) for UCDs is weak , 
dropp ing substantially after spectral class M7 ( Liebert et all Il999t Audard et al. . 
20071 ). Radio emission is expe cted to be und e tectable at the /iJy l evel, based 



on empirical X-ray scaling laws ( Guedel fc Bend . 1993 ; Benz fc Guedel . 1994 ). and 
the assumption that highly neutral UCD atmospheres are incapable of sustaining 
magnetic stresses which pervade the atmospheres of solar-type stars. However, a 
growing number of UCDs have been discovered that display significant emission 
at cm-wavelengths a nd suggest the presence of persistent kG-s c ale magnetic fields 



( Berger et all . boOlL iBurgasser fc PutmanL 120051 : lHallinan et all . 120081 : lOsten et al 



200' 



Furthermore, 
M0 



the quies c ent r a dio emission is nearly constant from spectral 



L5 ( Berger et all 2005: B erger . l2006h . It appears that for at least 



types M0 to 

some UCDs, the typical indicators of magnetic activity are not well correlated with 
decreased radio flux. 

Radio emission from these peculiar UCDs is typically broadband and unpolarized 
with high brightness temperature (10 8 — 10 9 K) during quiescence, and can 
10 11 



exceed 



Antonova et 



K with nearly 100% circular p olarization di 
fan. 120081 lHallinan et ai1 .l2008li. There is also 



olarization during bursts (jBergerl . [2002; 



evidence of long-term radio 
variability (jAntonova et all 120071). The a s sumed radiat i on me chanism was initially 



incoherent gyrosynchrotron (jBergerl l2002t I Berger et all 120051 ) from populations of 
mildly relativistic electrons with a power-law energy distribution. However, the 
high brightness temperature and high circular polarization seen during burst events 
sugges ted a coheren t radiation mechanis m such as the cyclo tron maser ins t abilit y 
(CMI. iMelrose et al|[l98llWingleel . ll985[ ). First suggested by lHallinan et al". I (l2006h . 
CMI has also been used to describe burst emi ssion from the polar, low density - high 



magnetic fiel d regio ns of magnetized planets (IZarkal 11998c lErgun et al 
dMutel et all Il998l ). and late- type flare stars (|Bingham et all l200ll : 



iity 

200ffl. Algol 



Kellett et al 



2002). The CMI model may also explain the quiescent emission, possibly created via 
depolarization of persistent maser sources (jHallinan et al. , 200G; Little fair et al.l , 2008; 



Yu et allliuTl . 



The exact mechanism(s) responsible for UCD radio emission (flaring on top of 
a quiescent background) is unclear, but it is likely due to a combination of both 
gyrosynchrotron emission and CMI. The emission frequency for each mechanism scales 
directly with the local magnetic field strength, emitting at the electron cyclotron 
frequency (Q C e,MHz — 2.8 • Bq). Gigahertz radio observations then require kG- 
scale fields, while detectable radio emission may exist at megahertz frequencies from 
regions of weaker field strength (~ 116 Gauss at 325 MHz). A vast majority of UCD 
observations have been limited to gigahertz frequencies where radio instruments are 
historically the most sensitive. Constraining the full spectral profile with addition of 
low frequency observations may provide valuable information to distinguish between 
the suggested emission processes, as well as reveal key differences in the atmospheres 
of UCDs compared to high mass M-dwarfs. Further, as UCDs are uniquely placed on 
the boundary between stars and sub-stellar objects, low frequency measurements also 
provide a guide for future observations of extrasolar planets. 
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2. Target History 



TVLM 513-46546 (hereafter TVLM513 ) is a spectral ty pe M9 dwarf with mass 
equaling 0.09 M sun and age > 1 Gyr (Reid et all . l2000h . Radio emissio n from 
this n ear-by source (~ 10.6 pc, Dahn et al.l ( 20021 )) was first detected by Berger 
(2002). Berger's observations at 8.5 GHz revealed both persistent stellar emission 
and occasional (2% — 10% duty cycle) pulses with flux densities exceeding 1 mJy. 
Detected bursts were highly circularly polarized and lasted multi ple minutes . Both 



the persistent and bu r st emission features we r e later confirmed bylOsten et al.l (12 006) , 
Hallirian et all (120061). lHallinan et al.l (|2007l ) . IBerger et all ((200a ). IForbrich fc Berger 
(120091). andlDovle et all (|2010|) throug h obse rvations at 8.5 GHz and 4.9 GHz. Further, 
Hallinan et all (|2007l ) and lLane et al.l (|2007l . optical) detec t ed a p ulse periodicity of ~ 
2 hr, consistent with the stellar rotation rate. lOsten et al.l (|2006l ) made a 5cr detection 
of TV1M513 at 1.4 GHz and i t is undete c ted a t lower frequencies at sensitivity levels 
similar to the GHz detections ( Antonova . 20071 ). 

2MASS J0036+1 821104 (2M0036 onwa rd) is a L3.5 br own dwarf with mas s 
- 0.06 - 0.074 M sun dSchweitzer et all. l200ll) age > 0.8 Gyr (|Burrows et all . l200ll ). 



and is l ocate d at 8.8 pc ( Dahn et all 120021 ). Radio observatons at 8.5 GHz by 



Bergerl (|2002) revealed quiescent emission with occasional bursts, si milar to those 
witnes sed on TVLM513. Subsequent measurements at 4.9 GHz by lHallinan et al.l 
( 20081 ) revealed strong (5x the quiescent level) circularly polarized pulses, lasting 5- 
20 min and with a 3 hr per iod corresponding t o the stellar rotation rate. However, 
observations at 8.5 GHz by IBerger et all (|2005l ); IForbrich fc Bergerl (|2009l ) detected 
no significant activity over multiple stellar rotation periods, indicating a possible 
absence of periodicity above 4.9 GHz. To date, there are no published observations of 
2M0036 below 4.9 GHz. 



3. Observations and Data Reduction 

Radio observations of UCDs TVLM513 and 2M0036 were conducted June 24 - 26, 
2007 using the NRAO Very large Array (VIA). Each source was observed for ~ 
10.5 hr (10 s integration) using 2 x 6.25 MHz bands centered at 327.5 MHz and 
321.6 MHz. Each frequency band was split into 15 spectral channels for the purpose 
of radio frequency interference removal and to mitigate bandwidth smearing. At 
observation time, the array utilized 23, 25 m antennas positioned in A configuration 
(maximum baseline ~ 35 km), resulting in a ~ 2.5 deg field of view and ~ 6" 
x 5" angular resolution. The radio flux density scale was set using amplitude 
calibrator 3C 286 (assumed 24.49 Jy at 327.5 MHz). Phase calibration and the 
receiver bandpass correction was performed using standard sources 1513+236 and 
0042+233 for TV1M513 and 2M0036 respectively. Measurements of each UCD and 
its corresponding phase calibrator were intertwined, performing 2 x 30 min primary 
target scans followed by a 5 min calibrator scan, then repeating. 

Data reduction and imaging were performed using both AIPS E and Obit @ 
software packages. The visibility data was calibrated using standard AIPS tasks. 
Automated data flagging, visibility self-calibration, and imaging was performed in 

I Astronomical Image Processing System, release 31DEC10 

§ Obit is developed and maintained by Bill Cotton at The National Radio Astronomy Observatory in 
Charlottesville, Virginia, USA and is made available under the GNU General Public License, version 
1.1.269-6-64b. 
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Obit. Lightcurves of TVLM513 and 2M0036 were made in both total (Stokes I) and 
circularly polarized (Stokes V) flux using AIPS task DFTPL. 

For comparison with published f.4 GHz quiescent emission observations of 
TVLM513, we analyzed VLA archival data from January 10, 2005 which contained 
an ~ 8 hr observation of 2M0036. The data consisted of 2 x 50 MHz bands centered 
on 1465 MHz and 1385 MHz obtained in the BnA hybrid configuration (maximum 
baseline ~ 21.2 km NS and 12.2 km EW). Amplitude and phase offsets were determined 
using calibrators 3C 147 (assumed 21.39 Jy at 1465 MHz) and 0042+233. Data editing, 
calibration and imaging were performed using standard AIPS routines. 



4. Results 



4.1. TVLM513 

We observe no significant unpolarized or circularly polarized radio emission associated 
with TVLM513 over the 10.5 hr observation. The p osition of T VLM513 is 
well determined (< 1"), based on closely spaced radio (jBerger et al. 



infrared (ICutri et al 



2008) and 



2003) detections and estimates of TVLM513's proper motion 



( Schmidt et all 120071) . Given the positional accuracy, we report a 2.5a non-detection 
limit to the total (Stokes I) and circularly polarized (Stokes V) quiescent radio flux. 
These limits are 795 fiJy for I and 743 /iJy for V. The relationship between our 2.5cr 
upper limit for TVLM513 and previously measured peak flux density values is shown 
in Figure Q] (scaled to 10.6 pc) and listed in Table [TJ 
Previous observations of TVLM513 indicate 



pulse period of ~ 1.96 hr 
( Hallinan et all 12006 ; Lane et al. . 2007 ). To search for burst emission and any 



potential periodicity, we constructed lightcurves at the known position of TVLM513 in 
both unpolarized and circularly polarized intensity with time resolutions between 10 s 
and 10.5 hr. We detect no significant variation to the measured radio intensity on 
any timescale. Also, no periodicity was found within the noise over ~ 5 stellar 
rotations, performing both a blind period search and by folding the flux values at 
the expected period. Lightcurves of the circularly polarized flux with 10 s and 5 min 
temporal resolution, along with the corresponding Lomb-Scargle periodogram of the 
10 s measurements are shown in Figure [2] Our 325 MHz non-detection implies a 
pulse flux density upper li mit of 2 4 mJy ( Stokes I) assum i ng a 1 0% duty cycle. High 
frequency observations by iBerger (|2002h : lHallinan et al] (j2007l ) observe pulse duty 
cycles in the range of 2-10%. We choose a 10% duty cycle to account for potential 
pulse dispersion at lower frequencies. 



4.2. 2M0036 

We searched the anticipated location of 2M0036 for polarized and unpolarized radio 
emission at observation frequencies of 325 MHz and 1.4 GHz. The position of 
2M0036 is known to an accuracy smaller than each of the synthesized beams (see 
Sec. [3J and is therefore well constrained. No significant radio emission was observed 
at either frequency. Our 325 MHz non-detection establishes a 2.5cx upper limit on 
the the quiescent flux < 942 /iJy in total intensity (Stokes I) and < 870 /^Jy in 
circularly polarized (Stokes V) intensity for the 10.5 hr observation. The non-detection 
at 1400 MHz (< 130 /iJy in Stokes I, < 95 /iJy Stokes V, 2. ho) sets an upper limit on 
the quiescent flux which is slightly lower than the predicted extrapolation from higher 
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Table 1. Measurement summary listing observing frequency v in GHz, 
observation length r in hr, and the recorded flux S in /xjy . Measurements are 
made using the VLA unless noted otherwise. 
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Burst 
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8.5 
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Burst 




23.09.2001 


8.5 


2 


135 ± 14 







frequencies when assuming a fiat spectrum (See Fig. |3|). 

No notable burst activity in the 325 MHz total flux or polarized flux measurements 
was observed on timescales from 10 s to 10.5 hr. The corresponding burst flux 
density upper limit is 2.8 mJy (Stokes I) assuming a 10% duty cycle We also 
detect no periodicity in the m easured flux oyer > 3 stellar rotations, assuming a 
stable rotation rate of 3.08 hr (|Hallinan et all l2008h . A lightcurve of the 325 MHz 
circular polarization measurements folded at the expected pulse period and a Lomb- 
Scargle periodogram of 10 s resolution data are shown in Figure HJ A similar flux 
variability search using the 1.4 GHz observation was not possible due to significant 
radio frequency interference which persisted for approximately 3 hrs during the session. 
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5. Discussion 

The absence of both flaring and quiescent emission from TVLM513 and 
2M0036 implies constraints on the plasma environments (i.e. electron energy, spacial 
distribution, etc.) of UCDs. Of the two sources, TVLM513 may offer the most 
stringent constraints, given the numerous GHz observations of stable pulse activity in 
the year preceding our June 2007 measurements (the closest burst is within 25 days, 



Doyle et all 12010). Here, we explore what limits can be inferred by our non-detections 



within the context of behavior seen on active stars and magnetized planets. 



5.1. Flare Emission Constraints 



Hallinan et al.l (|2008f ) have argued that the primary mechanism responsible for time 
variable emission from UCDs is the cyclotron maser instability (CMI). The CMI 
mechanism requires (1) a continuous source of energetic beamed electrons accelerated 
by a parallel electric field and (2) a converging magnetic field topology. As the electron 
beam propagates toward increasing field strength, conservation of the first adiabatic 
invariant causes transfer of parallel to perpendicular energy. The resulting velocity 
distribution f(v) becomes increasingly unstable (df/dv± > 0), leading to exponential 
wave growth at the electron cyclotron frequency. 

The CMI mechanism is quenched if the relativistic RX-mode cutoff frequency 
exceeds the cyclotron frequency. This constraint can be written in terms of an 
upper limit on the plasma j3, the ratio of electron plasma to cyclotron frequencies 
(|Mutel et al.Ll2006l) . 



P = 



< 



7-1 



Qce V 7 

with 7 representing the electron beam Lorentz factor. 



recast as 



/3 ~ 5 x 10 



-3 



< 1 



(1) 

For 7 — 1 <C 1, this can be 

(2) 



'G EkeV 

where n e is the electron density per cubic centimeter, Bq is the magnetic field in 
Gauss, and Ek e v is the beam energy in keV. 

By assuming a magnetic field configuration, an electron density profile, and a 
mean beam energy, this constraint can be mapped to a coronal volume for which 
the CMI mechanism is viable. We assume a dipole field with surface equatorial 
field stre ngth B a = 5 kG , cons i stent with estimates of 3 kG at the 8 GHz source 



altitude ([Reiners &, Basril 120071 : lHallinan et all 2008) but below the maximum of 
10 kG based on dynamo models ((BrowningF 2008 ) . The electron density profile is 
more problematic. Stellar coronal loop models suggest either a low-order power- 
law or expon ential dependence with radial distance, wit h a scale height of order 



the loop size ( Rosner et al. . 1978 : Collier Cameronl Il988f). We assume an example 



surface electron density n D = 10 cm 3 (jFludra et al. . 119991: lYu et~afl l201lh and an 



exponential dependence with a density scale height equal to the stellar radius. 

Figure [S£i illustrates the masing criteria (j3) in a model UDC magnetosphere for 
100 keV electron beams. Presumably, the periodic bursts observed on UCDs are 
generated in an active region extending along an co-rotating field line and emitting 
at a frequency consistent with the local magnetic field strength. Detections at both 
8.5 GHz and 4.9 GHz imply an active region extending to at least 1.5 stellar radii. This 
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is consistent with estimates from 2M0036 by lHallinan et al" (2008). Lower frequency 
observations probe larger radii, hence the highest observation frequency where (3 < 1 
and pulsed emission is absent defines the region's uppermost vertical extent. Assuming 
observations of 2M0036 and TVLM513 are typical of all UCDs and that the year-long 
history of pulse activity on TVLM513 continued throughout our 325 MHz observations 



(i.e. nearly 100% pulse duty cycle) , the absent pulse emission at 1.4 GHz (jOsten et al 



200G). 610 MHz (jAntonoval . 120071 ) and 325 MHz (this work, for strong flat spectrum 
bursts) suggests that the pulse source is confined within 2 stellar radii. 

Figure [SJd shows ratio of f3 to the critical f3 given in Equations 1-2 versus electron 
cyclotron frequency for a magnetic latitude of 65°. For 325 MHz, the CMI mechanism 
is only viable for electron beam energies E > 20 keV (E > 5keV at 1.4 GHz), whereas 
only above a few GHz beam energies below 1 keV can drive the CMI mechanism. 
This may also account for our non-detection: If CMI is driven by relatively low-energy 
beams, it could account for CMI-induced bursts at high frequencies, but suppression 
at lower frequencies. 



5.2. Quiescent Emission Constraints 

The multi-frequency observations (325 MHz - 8.5 GHz) are consistent with a flat 
gyrosynchrotron spectrum with indication in the 1.4 GHz measurements of 2M0036 of 
a possible low-frequency spectral break. While our 325 MHz observational limits are 
the most stringent to date, they are not sensitive enough to imply any real constraints 
on the quiescent emission mechanism for UCDs, requiring at least a factor of 5-10 
better sensitivity to rival the high frequency detections. The ability to make these 
needed measurements will rely on the completion of LOFAR (36 stations, ~ 300 /iJy 
2.5ct at 250 MHz in 10 hr) or the extension of the EVLA Low-Banc|J] retrofit to 
350 MHz (~ 100 ^Jy 2.5er in 10 hr with 15% fractional bandwidth). 



6. Summary and Conclusion 

We searched ultracool dwarfs TVLM 513-46546 and 2MASS J0036+1821104 for 
radio emission at 325 MHz. Measurements of the total (Stokes I) and circularly 
polarized (Stokes V) flux were made on timescales between 10 s and 10.5 hr. 
While strong continuous emission and multi-minute, circularly polarized pulses were 
previously recorded at 4.9 GHz and 8.5 GHz, no significant emission was measured at 
325 MHz. We set 2.5ct total flux limits of 795 /iJy and 942 /iJy for TVLM 513- 
46546 and 2MASS J0036+1821104 , respectively, for the quiescent emission over 
10.5 hr, consistent with a flat gyros ynchrotron spectrum with a potential spectral 



break near 2-3 GHz as indicated by ([Osten et al 



2006). Furthermore, we observe 



no significant variation in the 325 MHz radio flux at the expected location of each 
source. The absence of variable radio flux from TVLM 513-46546 below 1.4 GHz 
suggests that pulse emission from these UCDs (assuming a CMI electron acceleration 
model) originates from a source region confined below 2 stellar radii and/or is driven 
by electron beam energies less than a few keV. 

| The EVLA 74 MHz Low-Band system is scheduled to be completed Fall 2012 



Ultracool Dwarfs TVLM513 and 2M0036 



8 



7. Acknowledgments 

This paper utilizes data from VLA programs A0218 and AB1169. We would like 
to thank Bill Cotton for assistance with Obit. The National Radio Astronomy 
Observatory is a facility of the National Science Foundation operated under 
cooperative agreement by Associated Universities, Inc. This research was performed 
while the primary author held a National Research Council Research Associateship 
Award a the US Naval Research Laboratory. Basic research in radio astronomy at the 
Naval Research Laboratory is supported by 6.1 base funding. Part of this research was 
carried out at the Jet Propulsion Laboratory, California Institute of Technology, under 
a contract with the National Aeronautics and Space Administration. The LUNAR 
consortium is funded by the NASA Lunar Science Institute to investigate concepts for 
astrophysical observatories on the Moon. 



REFERENCES 



<) 



References 

Antonova, A. 2007, Ph.D. thesis, The Queen's University of Belfast 

Antonova, A., Doyle, J. G., Hallinan, G., Bourke, S., k Golden, A. 2008, A&A , 487, 
317. 10805 .45741 

Antonova, A., Doyle, J. G., Hallinan, G., Golden, A., k Koen, C. 2007, A&A , 472, 
257. IU7U7 . 06341 

Audard, M., Osten, R. A., Brown, A., Briggs, K. R., GMel, M., Hodges-Kluck, E., & 

Gizis, J. E. 2007, A&A , 471, L63. 10707 . 18821 
Benz, A. O., k Guedel, M. 1994, A&A , 285, 621 



Berger, E. 2002, ApJ , 572, 503. arXiv : astro-ph/0111317 



2006, ApJ , 648, 629. arXiv : astro-ph/0603176 



Berger, E., Ball, S., Becker, K. M., Clarke, M., Frail, D. A., Fukuda, T. A., Hoffman, 
I. M., Mellon, R., Momjian, E., Murphy, N. W., Teng, S. H., Woodruff, T., Zauderer, 



B. A., k Zavala, R. T. 2001, Nature , 410, 338. |arXiv: astro- ph/0102301 

Berger, E., Gizis, J. E., Giampapa, M. S., Rutledge, R. E., Liebcrt, J., Martin, E., 
Basri, G., Fleming, T. A., Johns-Krull, C. M., Phan-Bao, N., k Sherry, W. H. 2008, 
ApJ , 673, 1080. 10708 . TBTTI 

Berger, E., Rutledge, R. E., Reid, I. N., Bildsten, L., Gizis, J. E., Liebert, J., Martin, 
E., Basri, G., Jayawardhana, R., Brandeker, A., Fleming, T. A., Johns-Krull, C. M., 
Giampapa, M. S., Hawley, S. L., & Schmitt, J. H. M. M. 2005, ApJ , 627, 960. 
|arXiv : astro-ph/ 0502384 

Bingham, R., Cairns, R. A., k Kellett, B. J. 2001, A&A , 370, 1000 

Browning, M. K. 2008, ApJ , 676, 1262. 10712 ."16031 

Burgasser, A. J., & Putman, M. E. 2005, ApJ , 626, 486. |arXiv : astro-ph/0502365| 

Burrows, A., Hubbard, W. B., Lunine, J. I., k Liebert, J. 2001, Rev. Mod. Phys., 73, 
719 

Collier Cameron, A. 1988, M.N.R.A.S., 233, 235 

Cutri, R. M., Skrutskie, M. F., van Dyk, S., Beichman, C. A., Carpenter, J. M., 
Chester, T., Cambresy, L., Evans, T., Fowler, J., Gizis, J., Howard, E., Huchra, 
J., Jarrett, T., Kopan, E. L., Kirkpatrick, J. D., Light, R. M., Marsh, K. A., 
McCallon, H., Schneider, S., Stiening, R., Sykes, M., Weinberg, M., Wheaton, 
W. A., Wheelock, S., k Zacarias, N. 2003, 2MASS All Sky Catalog of point sources. 
(NASA/IPAC Infrared Science Archive) 

Dahn, C. C, Harris, H. C, Vrba, F. J., Guetter, H. H., Canzian, B., Henden, A. A., 
Levine, S. E., Luginbuhl, C. B., Monet, A. K. B., Monet, D. G., Pier, J. R., Stone, 
R. C, Walker, R. L., Burgasser, A. J., Gizis, J. E., Kirkpatrick, J. D., Liebert, J., 
& Reid, I. N. 2002, AJ , 124, 1170. |arXiv : astro-ph/0205050| 

Doyle, J. G., Antonova, A., Marsh, M. S., Hallinan, G., Yu, S., k Golden, A. 2010, 
A&A , 524, A15+ 

Ergun, R. E., Carlson, C. W., McFadden, J. P., Delory, G. T., Strangeway, R. J., & 
Pritchctt, P. L. 2000, ApJ , 538, 456 

Fludra, A., Del Zanna, G., Alexander, D., k Bromage, B. J. I. 1999, JGR , 104, 9709 

Forbrich, J., k Berger, E. 2009, ApJ Let. , 706, L205. IQ910". 13491 



REFERENCES 



10 



Guedel, M., & Benz, A. O. 1993, ApJ Let. , 405, L63 

Hallinan, G., Antonova, A., Doyle, J. G., Bourke, S., Brisken, W. F., k Golden, A. 
2006, ApJ , 653, 690. |arXiv : astro-ph/060 8556 

Hallinan, G., Antonova, A., Doyle, J. G., Bourke, S., Lane, C, k Golden, A. 2008, 
ApJ , 684, 644. 10805 .40101 

Hallinan, G., Bourke, S., Lane, C, Antonova, A., Zavala, R. T., Brisken, W. F., Boyle, 
R. P., Vrba, F. J., Doyle, J. G., k Golden, A. 2007, ApJ Let. , 663, L25. 10705 . 20541 

Kellett, B. J., Bingham, R., Cairns, R. A., k Tsikoudi, V. 2002, MNRAS , 329, 102 

Lane, C, Hallinan, G., Zavala, R. T., Butler, R. F., Boyle, R. P., Bourke, S., Antonova, 
A., Doyle, J. G., Vrba, F. J., k Golden, A. 2007, ApJ Let. , 668, L163. 10709 . 10451 

Liebert, J., Kirkpatrick, J. D., Reid, I. N., k Fisher, M. D. 1999, ApJ , 519, 345 

Littlefair, S. P., Dhillon, V. S., Marsh, T. R., Shahbaz, T., Martin, E. L., & 
Copperwheat, C. 2008, MNRAS , 391, L88. 10809 .21931 

Melrose, D. B., Dulk, G. A., & Hewitt, R. G. 1984, JGR , 89, 897 

Mutel, R. L., Menietti, J. D., Christopher, I. W., Gurnett, D. A., k Cook, 
J. M. 2006, Journal of Geophysical Research (Space Physics), 111, A10203. 
|arXiv : astro-ph/0 609802 

Mutel, R. L., Molnar, L. A., Waltman, E. B., k Ghigo, F. D. 1998, ApJ , 507, 371 

O sten, R. A., Hawley, S. L Bastian, T. S., k Reid, I. N. 2006, ApJ , 637, 518. 
arXiv : astro-ph/0509762 



Osten, R. A., Phan-Bao, N., Hawley, S. L., Reid, I. N., k Ojha, R. 2009, ApJ , 700, 
1750. 10905 .41971 

Reid, I. N., Kirkpatrick, J. D., Gizis, J. E., Dahn, C. C, Monet, D. G., Williams, R. J., 
Liebert, J., k Burgasser, A. J. 2000, AJ , 119, 369. |arXiv : astro-ph/9909336| 

Reiners, A., k Basri, G. 2007, ApJ , 656, 1121. |arXiv: astro-ph/0610365| 
Rosner, R., Tucker, W. H., k Vaiana, G. S. 1978, ApJ , 220, 643 
Schmidt, S. J., Cruz, K. L., Bongiorno, B. J., Liebert, J., k Reid, I. N. 2007, AJ , 
133, 2258. |arXiv : astro-ph/0701055| 

Schweitzer, A., Gizis, J. E., Hauschildt, P. H., Allard, F., k Reid, I. N. 2001, ApJ , 
555, 368. |arXiv : astro-ph/0103402| 

Winglee, R. M. 1985, JGR , 90, 9663 

Yu, S., Hallinan, G., Doyle, J. G., MacKinnon, A. L., Antonova, A., Kuznetsov, A., 
Golden, A., k Zhang, Z. H. 2011, A&A , 525, A39+. IT009. 15481 

Zarka, P. 1998, JGR , 103, 20159 



REFERENCES 



11 




10 13 l 1 1 I 

o.i i io 

Frequency (GHz) 

Figure 1. Peak luminosity of TVLM513 from all published radio observations, 
scaling the radio flux density to a distance of 10.6 pc. Diamonds indicate measured 
quiescent emission and stars indicate burst/pulse emission. Reported error bars 
are roughly the size of each marker. Bursts reported as lower limits are displayed 
with upward arrows. Non-detection limits (quiescent solid, 10% duty bursts 
dashed) are indicated with downward arrows ending at the la intensity. See 
Table [T] for specific radio measurements and author sources. 
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Figure 2. (top) Measured circularly polarized radio flux at the position of 
TVLM513 folded with a period of 1.96 hr. Dots indicate the 10 s resolution 
measurements. The solid line indicates the 5 min median value of the folded 
data, (bottom) Lomb-Scargle periodogram of the 10 s resolution flux values. 
Dashed lines indicate false alarm probabilities of 0.01 (99%), 0.1 (90%), and 0.5 
(61%). 
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Figure 3. Peak luminosity of 2M0036 from all published observations, scaling 
the measured radio flux density to a distance of 8.8 pc. Indicators are the same 
as in Figure [T] See Table [l] for specific radio measurements and author sources. 
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Figure 4. (top) Measured radio flux (Stokes V) at the position of 2M0036 folded 
with a period of 3.08 hr. Dots represent the 10 s resolution radio intensity. A 
solid line indicates the 5 min median value of the 10 s measurements, (bottom) A 
Lomb-Scargle periodogram of the 10 s resolution circularly polarized flux values. 
Dashed lines indicate false alarm probabilities of 0.01 (99%), 0.1 (90%), and 0.5 
(61%). 
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Figure 5. (a) Masing criteria above the surface of a ultracool dwarf, assuming 
a dipole field (B a = 5 kG), an exponentially decreasing electron density (n = 
10 9 cm -3 ), and an electron beam energy of 100 keV. Solid arcs are used to 
illustrate example magnetic field lines. Dashed curves indicate where the magnetic 
field magnitude is sufficient to produce cyclotron emission at 325 MHz, 610 MHz, 
1.4 GHz, 4.9 GHz, and 8.5 GHz. White to black shading indicates the value 
of beta (/3 = uj pe /Q C e)- (b) Plasma beta versus electron cyclotron frequency at 
magnetic latitude 65° for electron beam energies of 1 keV, 10 keV, and 100 keV. 
Beta has been normalized to its maximum allowed value for positive CM1 growth 
as given by equation (1). Note that at 325 MHz, only beam energies E > 10 keV 
result in CMI growth ((3 < 1), while above a few GHz, the CMI mechanism has 
positive growth at much lower energies. 



